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Interactions between glioma cells and their local environment are critical determinants of brain tumor growth, infiltration 
and neovascularisation. Communication with host cells and stroma via microvesicles represents one pathway by which 
tumors can modify their surroundings to achieve a tumor-permissive environment. Here we have taken an unbiased 
approach to identifying RNAs in glioma-derived microvesicles, and explored their potential to regulate gene expression in 
recipient cells. We find that glioma microvesicles are predominantly of exosomal origin and contain complex populations 
of coding and non-coding RNAs in proportions that are distinct from those in the cells from which they are derived. 
Microvesicles show a relative depletion in microRNA compared with their cells of origin, and are enriched in unusual or 
novel noncoding RNAs, most of which have no known function. Short-term exposure of brain microvascular endothelial 
cells to glioma microvesicles results in many gene expression changes in the endothelial cells, most of which cannot 
be explained by direct delivery of transcripts. Our data suggest that the scope of potential actions of tumor-derived 
microvesicles is much broader and more complex than previously supposed, and highlight a number of new classes of 
small RNA that remain to be characterized. 



Introduction 

Gliomas are the most common and deadly of adult primary 
brain tumors. They typically grow as widely infiltrative lesions, 
evade local innate immune surveillance, and induce prolifera- 
tion of local vasculature.' These features indicate that gliomas 
extensively modulate their local environment to achieve success- 
ful invasive growth. While tumor-stroma interactions are in part 
mediated by the release of soluble factors such as growth factors, 
there is growing interest in the role of small membrane-bound 
vesicles in mediating such interactions. 

Membrane vesicles are generally classified into three groups 
based on their size, surface membrane composition, and mecha- 
nism of release from the cell: exosomes, microparticles, and 
apoptotic bodies.^ Microparticles (100-1000 nm) and apop- 
totic bodies (1-5 p-m) are formed from cytoplasmic blebs that 
bud from the cell membrane and carry surface markers of their 
cell of origin. Exosomes (50—100 nm) carry distinct membrane 
proteins and are formed internally in multivesicular bodies that 
fuse with the plasma membrane to exit the cell. Microvesicles 
(MV), defined here as a mixed population of microparticles 



and exosomes, are capable of interacting with a variety of other 
cell types, and their target-cell specificity may be imparted by a 
distinct repertoire of membrane surface proteins.' ' Given their 
near ubiquity in multicellular organisms, MV appear to repre- 
sent a phylogenetically ancient form of intercellular communica- 
tion, relevant to normal biology as well as pathology (reviewed 
in refs. 6 and 7). 

Shedding of M V occurs in many (perhaps most) somatic cells 
in humans, and the process is increased in cancer cell lines and 
primary tumors.** Tumor-derived MV have been shown to sup- 
press the immune response, increase tumor progression, promote 
tumor invasiveness, and metastasis, and confer multidrug resis- 
tance.'" Gliomas produce MV that contain angiogenic proteins 
capable of stimulating endothelial cell growth and proliferation 
in vitro."" Glioma MV also carry mRNA, including that derived 
from mutant oncogenes, as well as microRNA and even DNA,'^ 
suggesting the scope and mechanisms of their actions on recipi- 
ent cells may extend beyond the promotion of local angiogen- 
esis. Furthermore, detection of glioma-derived MV in peripheral 
blood and cerebrospinal fluid suggest that they may be useful 
tumor biomarkers.""''^ 
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Figure 1. Characterization of microvesicles from U251 glioma cells. (A) 
Representative transmission electron micrograph of U251 microvesicle 
preparation. (B) Size distribution of a typical U251 microvesicle popula- 
tion as determined from electron micrographs. (C) Functional ontology 
analysis of microvesicle proteins identified by mass spectrometry. The 
proportion of microvesicle proteins in each ontology is shown on the left 
y-axis and the bar graph. Functional categories over-represented in the 
microvesicle proteins are presented in decreasing order significance with 
the lowest being P < 1 x 10 '.The right y-axis shows the statistical 
significance of each functional category (represented by the open 
diamonds). 



Small non-coding RNAs (sncRNAs) play a major role in 
the regulation of gene expression. The most abundant and best 
characterized group are microRNAs (miRNAs), which typically 
regulate the stability of mRNA transcripts and their translation 
into protein."* Other sncRNAs can regulate DNA transcription, 
chromatin structure, and genome integrity." SncRNAs have 
been identified in MV from various cell types, including glio- 
j^jg. 8,16,20,21 l^owever, most studies have focused on miRNA char- 
acterization, hence the full range and complexity of small RNA 
species in MV derived from gliomas is yet-to-be determined. In 
this study, we used microarray and deep sequencing to profile 
mRNAs and ncRNAs in glioma-derived MV. We find that spe- 
cific RNA species are selectively packaged in glioma MV, and 
that many of these are previously undescribed and of unknown 
function. Microvascular endothelial cells exposed to glioma 
MV exhibit significant gene expression changes, indicative of 



the capacity of these particles to modulate gene expression in 
recipient cells. 

Results 

Isolation and characterization of glioma microvesicles. 

Microvesicles (MV) purified from the media of the U251 glio- 
blastoma cell line are shown in Figure lA, as visualized by 
transmission electron microscopy (TEM). TEM revealed a rel- 
atively uniform population of membrane-bound vesicles rang- 
ing between 85—154 nm in diameter, with an average diameter 
of 103 nm (Fig. IB). Most MV exhibited a "cup" shaped 
morphology typical of exosomes under TEM.^^ MALDI-TOF 
mass spectrometry of the proteins extracted from MV prepara- 
tions identified a total of 112 proteins, of which 111 are known 
to be of exosomal origin (ExoCarta database,^^ exocarta.org; 
Table SI). Of the top 25 exosomal proteins in ExoCarta, 17 
were present in our MV preparation, including the tetraspanin 
CD9, an "exosomal marker" protein. Functional analysis of 
MV proteins using Ingenuity Pathway Analysis (IPA) identi- 
fied a large number of significantly enriched pathways, among 
which cancer- and apoptosis-related ontologies were prominent 
(Fig. IC). 

Overall, the size distribution, ultrastructural morphology, 
and protein composition strongly suggest that the glioma- 
derived MV in our preparations are largely exosomes. However, 
as we cannot discount the possibility that there may also be a 
contribution by other membrane-bound vesicles (such as small 
microparticles), we will continue to refer to the vesicles in this 
study using the collective term, microvesicles (MV). 

Glioma MV contain a complex population of small RNAs 
and an under-representation of miRNA. To characterize the 
small RNA content of glioma MV, we made small RNA librar- 
ies from RNA extracted from MV, as well as from the parent 
U251 glioma cells, and subjected them to SOLID sequenc- 
ing. Reads were mapped to the hgl9 version of the human 
genome using mapreads, annotated using custom scripts, and 
normalized for comparison (see Materials and Methods). The 
length distribution of small RNAs from U251 glioma cells 
(modal length 22 nt; Fig. 2A) was characteristic of the small 
RNA content of somatic cells and cancer cells, whose dominant 
small RNA species is miRNA,"' and the vast majority of these 
small RNA sequences were confirmed to be miRNAs (90.4%), 
with only very few reads mapping to other annotation catego- 
ries (Fig. 2B). In contrast, the size distribution of reads in MV 
was skewed toward smaller sequence length (modal length 20; 
Fig. 2A), suggesting an enrichment of other non-miRNA small 
RNAs. Annotation confirmed that a surprisingly modest propor- 
tion of MV small RNAs (38.7%) map to miRNAs, and there 
was a striking increase in the proportion of reads mapping to 
intergenic, repetitive, and intronic regions of the genome, as well 
as to tRNA and other miscellaneous RNA species (Fig. 2B). 
These findings indicate that glioma MV are depleted in miR- 
NAs relative to the cell from which they are derived, and that 
other, shorter small RNA species are selectively targeted to MV 
for export. 
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Figure 2. Specific miRNAs and miRNA variants are selectively packaged in glioma microvesicles. (A) Distribution of small RNA lengths in microvesicles 
and U251 cells. (B) Distribution of small RNA annotations in microvesicles and U251 cells. (C) Scatter plot of relative miRNA abundance in U251 parent 
cells and their microvesicles. Dashed lines indicate 2-fold enrichment/depletion; pinl< data points indicate enriched miRNAs with greater than 200 
counts in the microvesicles as well as greater than 2-fold enrichment in the microvesicles compared with the U251 cells; orange points indicate abun- 
dant miRNAs with greater than 5000 counts in the microvesicles; green indicates the overlap between these two groups. miRNAs that were present 
but did not meet the threshold of > 10 reads are represented with an abundance of 1. (D) Plots of three representative miRNAs that exhibit 5p/3p arm- 
switching in microvesicles relative to the U251 parent cells. 
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Table 1. Most abundant miRNAs in glioma microvesicles 



microRNA 


i^un 1 idii£cu ivi V 

reads 


ivuriiiciiizcu 1 

reads 


Fold enrich- 
ment in MV 


hsa-mir-21 


641159 


4352797 


0.15 


hsa-mir-99a 


28983 


88832 


0.33 


hsa-mir-23a 


27191 


95794 


0.28 


hsa-iTiir-30a 


20775 


45509 


0.46 


hsa-mir-30d 


20084 


63697 


0.32 


hsa-mir-30b 


19312 


41509 


0.47 


hsa-mir-22 


18736 


162910 


0.12 


hsa-mir-125a 


16350 


65083 


0.25 


hsa-let-7b 


16188 


53137 


0.30 


hsa-mir-25 


14731 


28718 


0.51 


hsa-iinir-221 


14722 


101775 


0.14 


hsa-iTiir-92b 


14560 


42586 


0.34 


hsa-mir-135b 


11213 


39642 


0.28 


hsa-mir-29a 


11084 


32056 


0.35 


hsa-mir-222 


11003 


25199 


0.44 


hsa-mir-100 


9704 


41055 


0.24 


hsa-mir-451a 


8131 


118 


68.91 


hsa-mir-4301 


8111 


636 


12.75 


hsa-mir-27b 


7868 


15314 


0.51 


hsa-mir-15b 


7836 


17002 


0.46 


hsa-mir-23b 


7620 


25071 


0.30 


hsa-mir-5096 


7613 


952 


8.00 


1 Ibd 1 1 1 1 1 JD/ O 


6575 


643 


10 23 


hsa-mir-30e 


6481 


16546 


0.39 


hsa-mir-374b 


6026 


33347 


0.18 


hsa-mir-339 


6004 


16376 


0.37 


hsa-mir-191 


5979 


42344 


0.14 


hsa-mir-4454 


5358 


2437 


2.20 



Specific miRNAs are enriched in MV relative to the par- 
ent glioma cells. Although miRNAs were relatively depleted 
in MV, a number of microRNA species were still abundant in 
both U251 glioma cells and their MV (Fig. 2C, Table 1), the 
most abundant being miR-21. This is unsurprising as miR- 
21 is one of the most consistently overexpressed microRNAs 
in cancer, including gliomas, and has been characterized as 
an "oncomir."^^ In addition, a small number of miRNAs were 
significantly enriched in the MV compared with parent U251 
cells, indicative of selective packaging. The most abundant of 
these are shown in Table 2A. Among those selectively pack- 
aged miRNAs, several have been reported in MV before: miR- 
451 was previously identified as a highly enriched miRNA in 
MV from HEK293 cells, endothelial cells and in breast cancer 
cell lines;^^'^** mir-1246 has also previously been reported to be 
selectively released in MV from breast cancer cell lines. Gene 
ontology analysis of the predicted targets of the top 10 selec- 
tively packaged miRNAs revealed significant enrichment of 
several functional pathways (Table 2B); most highly enriched 



was the KEGG pathway of glycosaminoglycan (GAG) biosyn- 
thesis, a pathway that plays a critical role in tumor growth and 
angiogenesis.^' 

In addition to the selective packaging of specific miRNAs, we 
also observed a number of miRNA processing variants in gli- 
oma MV compared with U251 cells (Fig. 2D; Table S2). These 
cases involved either a switch in arm bias when both 5p and 3p 
microRNAs were present, or a dominance of the "star" form of 
the miRNA (relative to the mature) in the MV. Considering that 
the seed sequences of 5p and 3p (and of mature and star) spe- 
cies will almost certainly be distinct, it is likely that the miRNA 
variants in the MV target a different repertoire of mRNAs to the 
miRNA in the glioma cells themselves. Differential packaging of 
miRNA variants derived from the same precursor also suggests 
that exosome biogenesis intersects at some point with miRNA 
biogenesis pathways. 

Transposon-derived small RNAs are prominent in glioma 
MV. In contrast to the U251 data set where a minority (5%) 
of small RNA reads mapped to repetitive elements, more than 
20% of reads in the MV data set mapped to repetitive elements 
(Fig. 2B), the majority of which were DNA transposons, or ret- 
rotransposons: LlNEs, LTR-type, and SlNEs (Fig. 3A). Increased 
transcriptional activity of transposable elements in cancer is well 
documented (reviewed in ref 30) and retrotransposon transcripts 
have previously been detected in microvesicles by microarray.^' 
We sought to determine whether these transposon-derived small 
(20 nt) RNAs exhibited signatures of processing (as opposed to 
being fragments of longer transcripts) by mapping to canoni- 
cal sequences representative of each retrotransposon class. With 
LlNEs we observed a fairly even distribution of reads mapping 
both sense and antisense with respect to the element, with the 
exception of a higher density around an antisense polyadenylation 
site known to be involved in gene breakage^^^ (Fig. 3B, arrow- 
head). The abundance of antisense LINE-derived small RNAs 
argues against these RNAs being a product of degradation. A 
canonical HERV-K sequence was used to determine the pattern 
of small RNA distribution for the LTR class of retrotransposons. 
Here we observed sense and antisense reads mapping along the 
open reading frames of the element, but the 5' and 3' LTR regions 
that act as promoters for these elements exhibited predominantly 
only antisense, or sense reads, respectively (Fig. 3C). Mapping to 
SINE elements (derived from 7SL sequence) produced the most 
striking pattern: the majority of reads mapped in a sense direc- 
tion to the distal 3' end of the element, and a smaller proportion 
concentrated in an antisense direction just downstream of the 
internal RNA Pol III promoter (Fig. 3D). Taken together, these 
data indicate that retrotransposon transcripts are processed into 
small RNA fragments for export out of the tumor cell. 

Novel small RNAs are selectively packaged in glioma MV. 
A large proportion of reads in the MV small RNA libraries were 
novel, and mapped to otherwise unannotated intergenic, as well 
as intronic locations (Fig. 2B). Loci with abundant reads showed 
"stacking" of sequences with clearly defined start and end posi- 
tions (representative examples are shown as insets in Fig. 4A— C). 
Their abundance and precisely defined distribution suggests that 
these novel intergenic and intronic reads are processed small 
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RNA species, as opposed to randomly distributed sequencing 
errors or degradation products. Furthermore, 50% (1056) of 
intron-derived small RNAs mapped exclusively in an antisense 
direction relative to the parent gene (e.g.. Fig. 4C), militating 
against these RNAs being a product of mRNA degradation or 
splicing. The majority of intergenic and intronic small RNAs 
were enriched in MV relative to the parent U251 cells (1976/2284 
(87%) and 1723/2092 (82%) with > 2-fold enrichment in MV 
relative to the parent U251 cells, respectively), indicating that 
most of these novel small RNAs are selectively packaged into 
MV. Interestingly, 1,904 intergenic and 1,668 intronic small 
RNA species were not detected in the parent glioma cells at all, 
occurring exclusively in the MV, some at very high abundance. 
This suggests that these specific small RNAs have functions only 
through MV trafficking. 

The functions of these intergenic and intron-derived small 
RNAs are at present unknown, but they are unlikely to be novel 
microRNAs: they have an average length of 20 nt, and we did 
not detect any significant hairpin structures typical of miRNA 
precursors in the surrounding sequence. 

Microvesicles are highly enriched in vault RNAs. Some of 
the most abundant and enriched small RNA species in our MV 
data set were 20mers that mapped to loci encoding an RNA 
family known as "vault" RNA; these small RNAs were virtu- 
ally absent from the parent U251 cells (Fig. 5A). Vault RNAs 
complex with three proteins to form the vault organelle, and 
although the function of these RNAs is not understood, the vault 
organelle has been implicated in multidrug resistance in cancer.'^ 
Vault RNA has recently been shown capable of being processed 
by Dicer into six distinct small vault RNAs (svRNAs),^'' and 
at least some of these have been shown to exhibit miRNA-like 
Ago2-dependent repression of semi-complementary targets. In 
our data set we observed svRNAs mapping predominantly to 
the 3' end of the vauk transcripts VTRNAl-1 and VTRNA2-1, 
although for VTRNAl-1 (the most abundant) we also observed 
reads mapping from the 5' end (Fig. 5B). This contrasts with the 
original description of svRNAs, in which the majority of svR- 
NAs were derived from the 5' end.'^ We did not detect the pro- 
tein components of the vault complex in our mass spectrometry 
analysis of MV proteins, but the coverage of small RNA reads 
at VTRNAl-1 prompted us to perform northern blots to assess 
whether full-length vault RNAs were also present in the MV. We 
found that full-length VTRNAl-1 was very abundant in total 
RNA from the MV, much more so than in the total RNA from 
the parent U251 glioma cells (Fig. 5C). We did not observe the 
small vault RNA species on the northern blot, indicating that 
full-length vault RNA is present in MV in vast excess to its pro- 
cessed forms. 

The presence of full-length vault RNA in MV prompted us to 
then ask what other longer RNAs are present in glioma MV. We 
compared the expression profiles of MV and their parent U251 
cell cultures by hybridizing total RNA from triplicate MV/U251 
preparations to Affymetrix 1.0 ST arrays. Overall, there was a 
good correlation between the RNA profiles of MV and their par- 
ent U251 cells (r- = 0.88; Fig. 6A). However, a subset of tran- 
scripts was selectively enriched in MV; those most enriched (> 



Table 2A. Selectively packaged mlRNAs in glioma microvesicles" 



microRNA 


FnlH Anrirhmpnt 

ruiu CI II 1^1 ■ ■ ■ ici 1 L 

inMV 


NnfiTiAliTpH 

MV reads 


NnfiTiAliypri 

IHwl IIICII l£VU 

U251 reads 


hsa-mir-451a 


69 


8131 


118 


hsa-mir-4301 


13 


8111 


636 


hsa-mir-5096 


8 


7613 


952 


hsa-mir-3676-5p 


10 


6575 


643 i 


hsa-mir-4454 


2 


5358 


2437 


hsa-mir-1303 


73 


2350 


32 


hsa-mir-1273a 


78 


2100 


27 


hsa-mir-619 


5 


1954 


398 i 


hsa-mir-448 


3 


1799 


520 


hsa-mir-1246 


9 


914 


98 


hsa-mir-4792 


666 


666 


< 10 


hsa-mir-5095 


9 


615 


1 


hsa-mir-1273 g 


610 


610 


< 10 


hsa-mir-4256 


2 


502 


218 


hsa-mir-4255 


11 


376 


34 


hsa-mir-5100 


9 


313 


34 


AC068946.1 


270 


270 


< 10 


hsa-mir-1285-1 


11 


266 


25 


hsa-mir-1269b 


248 


248 


< 10 


hsa-mir-4500 


227 


227 


< 10 


hsa-mir-1273d 


207 


207 


< 10 


hsa-mir-4443 


171 


171 


<10 



"With > 2-fold enrichment in iVlV, as shown in Figure 2C. 



Table 2B. KEGG pathways targeted by the top 10 selectively packaged 
miRNAs 



KEGG pathway 


P value** 


Glycosaminoglycan biosynthesis - heparan sulfate 


5.077E-26 


Endocrine and other factor-regulated calcium reabsorption 


2.920E-03 


ECM-receptor interaction 


9.994E-03 


''DIANA-Mirpath, with Benjamini-Hochberg correction. 





3-fold) are shown in Figure 6B. Consistent with our observations 
with northern blotting (Fig. 5C), we again find that VTRNAl-1 
is the most enriched transcript in MV. Gene ontology analysis of 
these MV transcripts obtained a number of enriched functions 
related to cell morphology, growth, signaling, and development 
(Fig. SI). 

Glioma MV induce gene expression changes in cultured 
vascular endothelial cells. We exposed cultured human brain 
microvascular endothelial cells (EC) to glioma MV to determine 
whether the molecular cargo of MV could alter the transcriptome 
of recipient cells. After 24 h of exposure the EC were washed 
and RNA harvested for expression profiling using Affymetrix 
1.0 ST arrays; RNA from untreated EC of the same confluence 
were used for comparison and the experiment was performed 
in triplicate. We found 54 transcripts that had a greater than 
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Figure 3. Repeat element-derived small RNAs are abundant In microvesicles. (A) Abundance of repeat-derived small RNAs in microvesicles and U251 
cells. The distribution of small RNA reads across a canonical element (sense, above; antlsense, below) is shown for (B) LINEs, (C) LTRs, and (D) SINEs. 



2-fold expression change in EC upon exposure to MV: 35 were 
upregulated, wtiile 19 were downregulated. Of the 35 upregu- 
lated transcripts only four corresponded to transcripts that were 
abundant in MV (Fig. 6C; Fig. S2), indicating that the increases 
in EC gene expression are unlikely to be due to direct delivery 
of mRNA from the MV in most cases. Of the genes that were 
downregulated in the EC in response to MV exposure, many are 
predicted targets of microRNAs that are abundant or selectively 
packaged in the MV (Fig. 6D). 

Interestingly, we also observed that the majority (17/20) of 
genes that harbored antlsense intron-derived sRNAs in the 
MV (Fig. 4C) were also downregulated in EC exposed to MV 
(Table 3), suggesting that these intronic antlsense small RNAs 
may be able to modulate their parent gene's expression when 
delivered to recipient cells. 

Discussion 

In this study, we have found that microvesicles (MV) released 
from U251 glioma cells are predominantly of exosomal origin 
and carry a complex cargo of RNA species, many of which are 
novel and of unknown function. Taking an unbiased approach 
with small RNA deep sequencing, we have found that the reper- 
toire of small noncoding RNAs in MV is strikingly distinct from 



the repertoire in their cell of origin, indicating that many RNA 
species are produced within tumor cells specifically for export 
out of the cell. Exposure of brain microvascular endothelial cells 
to glioma MV resulted in significant changes in endothelial cell 
gene expression, which reflects the capacity of these particles to 
mediate intercellular communication. Our data suggest that glio- 
mas can modify the transcriptional landscape of their local envi- 
ronment through a variety of microvesicle-delivered RNA-based 
pathways. Many of these RNA-based mechanisms remain to be 
characterized. 

Most previous studies on the RNA content of MV have 
focused on miRNAs;'^'"''^''" these small (-22 nt) non-coding 
RNAs are post-transcriptional regulators formed by the process- 
ing of a larger precursor transcript with a hairpin intermediate.'* 
With our unbiased approach, we have found that miRNAs are 
relatively depleted in MV, although a small subset is specifically 
enriched (Table 2B). The most enriched miRNA, miR-451, 
plays an important role in glioma cell proliferation and migra- 
tion,'"' but it has been reported to be selectively enriched in MV 
from non-neoplastic cell lines as well,^^'*' indicating that selective 
export of this particular miRNA is not tumor-specific. However, 
available evidence suggests that selective export of most other 
miRNA species in glioma MV is unique; this points to potential 
tumor-specific modulatory roles of these miRNAs. 
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Figure 4. Novel small RNAs derived from Intergenic 
regions and Introns are selectively packaged in glioma 
microvesicles. (A) Bar graph showing the abundance of the 
top 20 intergenic loci producing small RNAs found almost 
exclusively in microvesicles. Inset histogram shows the 
characteristic stacked distribution of reads at a representa- 
tive locus. (B) Bar graph showing the abundance of small 
RNAs mapping to introns in a sense direction to the gene 
shown on the y-axis. Inset histogram shows the characteris- 
tic stacked distribution of reads at representative locus. (C) 
Same as in (B) except the small RNAs map in an antisense 
direction to the genes indicated on the x-axis. Note that 
these antisense intronic small RNAs are found almost exclu- 
sively in the microvesicles and not the parent U251 cells. 
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Most miRNAs show predominant accumulation of 
one "arm" of the precursor transcript after processing, 
and this is generally thought to be the mature functional 
effector miRNA, stabilized by incorporation into the 
RISC complex.'^ In a number of cases we found that the 
dominant/mature miRNA in MV was derived from the 
opposite arm of the miRNA precursor to that which was 
dominant in the parent U251 cells. This surprising find- 
ing suggests that miRNA biogenesis and endosomal/exo- 
somal processing may be interlinked, and that both arms 
of miRNA precursors are utilized more commonly than is 
currently appreciated, albeit in different contexts. 

A previous study on glioma MV used antibody arrays 
to identify angiogenic proteins (such as angiogenin, IL-6 
and IL-8) within MV that were capable of stimulating 
endothelial cell growth and tubule formation in vitro."' 
We did not find these angiogenic proteins in our mass 
spectrometry analysis of glioma MV proteins, but we did 
find that many transcripts downregulated in endothe- 
lial cells exposed to glioma MV are predicted targets of 
miRNAs that are either abundant or specifically enriched 
in the MV. KEGG pathway analysis of these predicted 
miRNA targets identified GAG synthesis, Ca^* reabsorp- 
tion, and ECM-receptor interactions as highly significant, 
over-represented functions. All of these pathways are rel- 
evant to glioma biology: GAG biosynthesis is a power- 
ful modulator of tumor growth and angiogenesis,^' Ca^* 
signaling modulates cell death and proliferation,*^ and 
ECM-receptor interactions are critical for tumor infiltra- 
tion.*^ Our findings suggest that modulation of vascular 
endothelial cell behavior by MV can also occur through 
small RNA pathways, in addition to protein-based 
mechanisms. In addition to any sequence-specific effects 
of MV miRNA, recent reports suggests that these (and 
perhaps other) small RNA species may activate Toll-like 
receptors in a sequence-independent manner, leading to 
multiple downstream immunomodulatory events in the 
recipient cells.*''''''' 

Perhaps the most striking finding is the abundance 
of novel small RNAs in glioma MV. Nearly half of the 
small RNA population identified are novel small RNAs 
that map to otherwise unannotated intronic and inter- 
genic regions. The intergenic RNAs are particularly intriguing: 
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Figure 5. Glioma microvesicles are highly enriched in vault RNAs. (A) Bar graph 
showing abundance of small RNAs derived from vault RNA family members 
in microvesicles; these processed vault RNAs are absent from the parent U251 
cells. (B) Histograms showing the distribution of microvesicle small RNAs across 
VTRNA1-1 and VTRNA2-1.(C) Northern blot of total RNA from microvesicles and 
U251 cells using a probe against the 3' end of VTRNA1-1. Mature vault transcript 
is indicated by the arrow. 



many are found in gene-poor regions whose sequence conserva- 
tion extends only to primates. Tire function of these RNAs is 
completely unknown. Intronic-derived sRNAs are found both 
sense and antisense to the parent gene, and exhibit processing 
features that indicate they are not simply degradation products 
of their parent transcripts. Furthermore, the overwhelming 
majority are undetectable in the cells from which the MV are 
derived, suggesting that the primary function of these small 
RNAs is achieved via MV-trafficking. In the case of antisense 
intron-derived reads, some function can be inferred by the 
endothelial cell co-culture experiments, where parent tran- 
scripts in the recipient cells are downregulated (Table 3); this 
suggests a role for these small RNAs in gene repression that is 
not intended in the cell in which the small RNA is produced. 



Repeat-derived small RNAs are the second most abun- 
dant class of small RNA in glioma MV and at least in the 
case of SINEs and LTR elements, the processing features 
of these small RNAs suggests that they are not the prod- 
uct of degradation. Retrotransposon RNAs have previ- 
ously been described in glioma MV, although they were 
detected by microarray and presumed to be full-length 
transcripts.^' While it is possible that the small RNA frag- 
ments may have contributed to the microarray signals in 
that study, it is also possible that, like the vault RNAs, 
both full-length and processed RNAs are packaged into 
glioma MV. The function of these small RNAs is again 
completely unknown; however, similar retrotransposon- 
derived small RNAs have been observed in murine oocytes 
and pre-implantation embryos,''^'" where there is some 
evidence that they are involved in gene silencing during 
early development.^^ It is tempting to speculate that these 
RNAs may have silencing roles in host cells exposed to 
MV that could promote a tumor-permissive environment; 
however, given the vast numbers of retrotransposons 
interspersed among all the genes in the genome, it may be 
impossible to identify specific targets. 

Vault complexes are very large (13 MDa) ribonucleo- 
protein organelles, three times the size of ribosomes, 
whose function is not well characterized.'"* They are com- 
posed of three proteins and several non-coding RNAs. 
The major vault protein accounts for 70% of the mass of 
vaults, and has been implicated in multidrug resistance 
in a variety of tumors, including gliomas.*^ Although it 
is abundant in glioma cell lines, we did not identify the 
major vault protein (or other vault proteins) in glioma 
MV by mass spectrometry. We did however identify 
an abundance of full-length vault RNA, VTRNAl-1, 
in MV by microarray and northern blotting, and deep 
sequencing revealed a large proportion of small RNA 
reads that were consistent with processed vtRNA. Both 
full-length and small vtRNAs were highly enriched in 
MV relative to the U251 parent cells; similar processed 
VtRNAs have also been identified in vesicles obtained 
from the media of murine T-cell and dendritic cell co- 
cultures.'" The processing of vtRNA into functional 
fragments has been demonstrated in breast cancer cells, 
but in contrast to this report, we did not see repression of pre- 
dicted targets of small vtRNAs in the EC exposed to the MV. 
Thus, it is likely that small vtRNAs have additional targets or 
different functions yet to be characterized. Taken together with 
the abundance of full-length vtRNA in the absence of major 
vault protein, vault RNAs and their derivatives likely have func- 
tions in MV separate to the vault ribonucleoprotein complex. 

The study of tumor MV is relatively young, but those 
studies that have focused on their RNA content have already 
shown that MV-derived miRNAs are likely to play a signifi- 
cant role in regulation of gene expression in host cells (reviewed 
in refs. 39 and 51). The most striking finding of this study is 
that glioma MV contents are relatively depleted in miRNAs 
and instead are highly enriched for unusual or completely novel 
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Figure 6. Glioma microvesicles alter the transcriptional profile of exposed endothelial cells. (A) Scatterplot comparing transcriptional profiles of U251 
cells and their microvesicles. (B) Heat map showing transcripts that are enriched in microvesicles relative to U251 cells (> 1 .5 log^ enrichment, FDR < 
0.15). (C) Heat map showing the relative expression levels of transcripts that are significantly upregulated in endothelial cells in response to exposure to 
glioma microvesicles. Transcripts which are in common between B and C are colored. (D) Heat map showing transcripts that are significantly downregu- 
lated in endothelial cells in response to microvesicle exposure. Transcripts that are predicted targets of miRNAs in the microvesicles are shown in red. 



non-coding RNAs, including small RNAs from primate-spe- 
cific regions of the genome. The functions of most of these are 
yet to be characterized, but point to unexpected complexities 



in intercellular communication mediated by membrane-bound 
vesicles. Microvesicle RNAs represent a new and intriguing 
mechanism of tumor-host interactions that may provide novel 
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Table 3. Response of EC genes to MV exposure where the gene has 
antisense intron-derived sRNAs in MV 



Gene with MV-associated 
antisense sRNA 


Log^fold change in EC+IVIV 


LRPAP1 


-0.012 


C4orf37 


-0.054 


TMC03 


-0.190 


ST8SIA1 


-0.100 


THSD7B 


-0.146 


G6PC2 


-0.272 


VIT 


-0.252 


IL1RAPL2 


0.008 
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-0.056 


RANBP3 


0.016 


LRP1B 


-0.004 


PDE11A 


-0.088 


WDR64 
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ABCA13 


-0.165 


RAD18 


-0.094 


EPHA3 


-0.086 


GALNTL6 


-0.102 


PIGG 


-0.139 


RBMXL3 


0.162 


TMPRSSllF 


-0.154 



therapeutic opportunities in the treatment of ghoma and other 
tumors. 

Materials and Methods 

Cell culture and MV preparation. Human Glioblastoma 
(Astrocytoma) cells (U251-MG; U251) were obtained from the 
ATCC. They were cultured in RPMI-Glutamax (Invitrogen) with 
10% fetal calf serum, 1% glucose, and 1% penicillin/streptomycin 
filtered through a 0.1 |xM filter. U251 cells were seeded at 2 x lO*" 
cells in a T75 culture flask and media collected from cultures after 
the cells reached 80-90% confluence. Media was centrifuged at 
1800 g for 10 min to remove cellular debris. The supernatant was 
then collected and centrifuged at 18000 g for 45 min at 15 °C. 
Pelleted microvesicles were resuspended in 500 |jlL of phosphate 
buffered saline, centrifuged again at 18000 g for 45 min and the 
pellet snap-frozen in liquid nitrogen. Microvesicle samples were 
verified by flow cytometry analysis, as described previously,'^ and 
electron microscopy as described below. Multiple microvesicle 
preparations were then pooled for RNA and protein extraction. 
From a total of 225 mL of U251 conditioned media we obtained 
approximately 1.3 x lO"" MV; this number of MV yielded a total 
of 2.36 |xg RNA. The contribution of any bovine MV from the 
fetal calf serum remaining after filtration was determined to be 
minimal: identical preparations from equivalent amounts of media 
alone yielded no obvious MV pellet, and less than 50 ng of RNA. 

Transmission electron microscopy. Microvesicles were har- 
vested as above, and the final pellet was resuspended in 2.5% 



glutaraldehyde in 0.1 M cacodylate buffer pH 7.4 and fixed for 
24 h at 4°C. The sample was centrifuged at 3000 rpm for 20 min 
and the pellet was processed routinely for electron microscopy. 
Ultrathin sections were stained with uranyl acetate and lead 
citrate and examined on a Philips CMIO electron microscope. 

Protein analysis. Microvesicle pellets were denatured in LDS 
Sample Buffer (Life Technologies) diluted in PBS and 100 mM 
DTT Samples were heated to 95 °C for 5 min and centrifuged at 
17000 g to precipitate any insoluble material and separated using 
a NuPAGE® Novex 4-12% Bis-Tris Gel (Life Technologies). The 
entire lane was then excised from the gel and sent to the Bioanaly tical 
Mass Spectrometry Facility at the University of New South Wales 
for protein identification. Scaffold software (Proteome Software'^) 
was used to probabilistically validate protein identifications from 
Mascot scores using the XITandem and ProteinProphet computer 
algorithms. Validated proteins were defined as those with at least 
five peptides detected with a minimum identification probability of 
95%, and a minimum protein identification probability of 99.9%. 
Functional ontology analysis of validated microvesicle proteins was 
performed using Ingenuity Pathways Analysis software (Ingenuity 
Systems; www.ingenuity.com). Right-tailed Fisher's exact test was 
used to calculate a p value determining the probability that each 
biological function had occurred due to chance alone, followed by 
Benjamini-Hochberg correction for multiple testing. 

Small RNA sequencing. Total RNA was extracted from 
frozen microvesicle and U251 cell pellets with TRIzoF (Life 
Technologies) according to the manufacturer's instructions. 
Small RNA libraries were generated using the NEBNext® Small 
RNA Sample Prep 3 kit (New England Biolabs) according to kit 
protocols. Libraries were sequenced on an Applied Biosystems 
SOLID 5500 X 1 using 35 bp chemistry. 

Small RNA mapping. SOLID color space reads (35mers) 
were mapped to the human genome using LifeScope (Life 
Technologies). First reads were filtered through a default filter 
file using a 25 base seed allowing up to three errors (25.3.0) iden- 
tifying tRNA, rRNA, and repetitive elements, including LINE, 
SINE, and LTR elements. Tags were then mapped to mirBase 
vl8 (18.2.0). Remaining reads were mapped to human genome 
build 19 (20.1.0). 

Annotation of mapped reads. Mapped reads were annotated 
using custom Perl scripts and the following databases: mirBase 
vl8, UCSC hgl9 Repeat Masker track, UCSC piRNA anno- 
tations and the Ensembl database version 67. Each read was 
assigned a gene accession, a biotype/RNA class, and orienta- 
tion if its coordinates were found to overlap with those recorded 
in one of the databases. This was performed in a hierarchical 
manner where the mirBase data set was queried first, followed 
by the Repeat Masker track, the piRNA track, and finally the 
Ensembl database through a local installation of the ensemble 
API. Reads with no known annotations in these data sets were 
assigned as "intergenic." Reads mapping to or near protein cod- 
ing genes were further grouped into the following subclasses 
upstream/downstream, exonic, intronic, antisense, 5' UTR, or 
3' UTR. 

Counts and enrichment analysis. All counts were expressed 
as normalized values, per 10 million mapped reads, prior to 



1342 



RNA Biology 



Volume 10 Issue 8 



filtering. A two-pass filter was applied to each of the reads befiare 
being included in further analysis. First, reads were required to 
have a Map Quality (MAPQ) of 10 or more before progressing 
to the second filter step where the number of reads starting at a 
given map position should be greater than or equal to 10. After 
filtering the number of mapped reads was 2,856,057 for U251 
cells and 1,725,859 for MV. 

Repeat element analysis. Reads mapping to three of the major 
classes of repeats LINE, SINE and LTR were remapped against 
a single representative full-length sequence of that repeat element 
with SHRIMP^'' using a 14 base seed and allowing up to five mis- 
matches to account for their characteristic degeneracy. The elements 
chosen for remapping were LINE L19092, SINE7SL NR002715 
and LTR HERV K JN675068. The resuhing alignments were 
plotted against a schematic of the representative sequence, separat- 
ing tags mapping to the positive or negative strand. 

Co-cultures. Primary human brain microvascular endothelial 
cells (EC) were obtained from Cell Systems, Inc. EC were cul- 
tured on plasticware coat with 3% collagen and were maintained 
in EBM-2 complete media (Lonza) with 5% fetal calf serum, 
ascorbic acid (5 JJtg/ml), hydrocortisone (Sigma, 1.4 |JLmol/L), 
chemically defined lipid concentrate (Invitrogen, 1:100 dilution), 
HEPES (Invitrogen, 10 mmol/L) and (3-Fibroblast Growth Factor 
(Sigma, 1 ng/ml). EC were used from four to passage 11. Purified 
U251 MV were added to EC cultures at -50% confluence at a 
ratio of lOMV/EC. EC were co-cultured with MV for 24 h and 
then washed extensively with PBS prior to harvesting for RNA 
extraction. 

Microarray. Total RNA was extracted from frozen microves- 
icles, U251 cell or endothelial cell pellets with TRIzol® (Life 
Technologies) according to the manufacturer's instructions 
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